Collective flow in p-Pb and d-Pb collisions at TeV energies 
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We apply the hydrodynamic model to the dynamics of matter created in p-Pb collisions at y^sjvjv = 
4.4TeV and d-Pb collisions at ^/sNN ~ S.llTeV. The fluctuating initial conditions are calculated in 
the Glauber Monte-Carlo model for several centrality classes. The expansion is performed event by 
event in 3 + 1-dimensional viscous hydrodynamics. Noticeable elliptic and triangular flows appear 
in the distributions of produced particles. 



(N 

T— I 

O 
(N 

a. 

oo: 

(N- 

CD. 

fi: 
>: 

in. 
^ . 

On- 

o: 



PACS numbers: 25.40.Ep, 25.45.De, 13.85Hd, 25.75.Ld 

Keywords: hadron-nuclcus collisions, viscous hydrodynamic model, collective flow 



X: 



I. INTRODUCTION 

The large multiplicity of particles emitted from the 
small interaction region in relativistic heavy-ion collisions 
implies that a fireball of very dense matter is formed. 
Experiments at the BNL Relativistic Heavy Ion Collider 
(RHIC) and the CERN Large Hadron CoUider (LHC) 
[il, |2|| have demonstrated the appearance of a collective 
flow in the expanding fireball. The physical picture is ex- 
pected to be different in the interaction of small systems, 
proton-proton, proton-nucleus or deuteron-nucleus. At 
RHIC energies the density of matter created in d-Au in- 
teractions is small and does not cause jet quenching [l| . 
d-Au and p-p interactions are treated as a baseline ref- 
erence to evidence new effects in nucleus-nucleus colli- 
sions, beyond a simple superposition of nucleon-nucleon 
(NN) collisions. With the advent of proton-proton colli- 
sions at several TeV center of mass (cm.) energies at the 
LHC, it has been suggested that some degree of collective 
expansion appears in high multiplicity p-p events 
However, no direct experimental evidence exist for such 
a collective expansion in p-p interactions. 

At the LHC, p-Pb collisions can be studied in the fu- 
ture, experiments with d-Pb or other asymmetric systems 
are also possible, but with additional technical difficulties 
Q . Estimates of the hadron production in p-Pb interac- 
tions at TeV energies take into account nuclear effects 
on the parton distribution functions, saturation effects, 
but do not assume the formation of a hot medium 0, Q . 
Experiments with p-Pb beams should provide an input 
for models used in heavy-ion collisions for the calculation 
of dense medium effects on hard-probes. 

The expected multiplicity and size of the interaction 
region in central p-Pb and d-Pb collisions at TeV en- 
ergies are similar as in peripheral (60 — 80% centrality) 
Pb-Pb collisions at ^/snn = 2.76TeV 9]. This raises 
the possibility that hot and dense matter is formed in 
such collisions. For strongly interacting matter, the as- 
sumption of local equilibrium is a good approximation 
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and relativistic hydrodynamics can be used to follow the 
evolution of the system (loj . Quantitative predictions for 
the elliptic flow have to account for finite deviations from 
local equilibrium in the rapidly expanding fluid [lll - fisj . 

In order to test the assumption of the formation of a 
dense fluid in p-Pb and d-Pb interactions and to esti- 
mate possible effects of its collective expansion, we apply 
the viscous hydrodynamic model to calculate the spectra 
of emitted particles. The goal of this study is to have a 
quantitative prediction of the elliptic and triangular flows 
and of the transverse momentum spectra for comparison 
with future experiments. The dynamically evolved den- 
sity of the fireball from hydrodynamic simulations can be 
used in the calculations of the parton energy loss in such 
small systems. 

The task requires the use of the most sophisticated ver- 
sion of the hydrodynamical model: event by event 3 -f 1- 
dimensional (3 -f 1-D) viscous hydrodynamics. While a 
good description of many collective phenomena in heavy- 
ion collisions can be obtained in the perfect fluid hydro- 
dynamics m2 + 1-D [l3,[ll| or 3 -M-D [H], to calculate 
the azimuthally asymmetric flow in small systems such 
as p-Pb or d-Pb collisions one has to use viscous hy- 
drodynamics. In collisions of symmetric nuclei 2 -I- 1-D 
boost-invariant viscous hydrodynamics is routinely be- 
ing applied for observables at central rapidities [111, ■ 
In p-Pb or d-Pb interactions the energy density and the 
final particle distributions depend strongly on rapidity. 
This forces the use of 3 -f 1-D hydrodynamics to obtain 
realistic particles spectra at different rapidities. Only 
recently 3 -I- 1-D viscous hydrodynamic simulations be- 
came available [l^ [l^. In proton or deuteron interac- 
tion with a nucleus the shape of the interaction region 
fiuctuates widely from event to event. Unlike in inter- 
actions of heavy ions, using the average density is not a 
reliable approximation. Event by event 3 -I- l-D perfect 
fluid hydrodynamics is used by several groups jl7| . The 
inclusion of event by event fluctuations is important in 
the description of the initial eccentricity and triangularly 
of the fireball [l^, [T7l - [2l| . Only one group is using an 
event by event 3 -I- 1-D viscous hydrodynamic code for 
heavy-ion collisions [H, HI] . 

As the size and the life-time of the system decrease the 
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hydrodynamic model becomes less justified. A sizable el- 
liptic flow is observed in peripheral Pb-Pb collisions at 
the LHC, which proves that substantial rescattering oc- 
curs in the evolution of the fireball. By itself it does not 
prove that the hydrodynamic regime is applicable in such 
collisions, as some elliptic fiow can be generated through 
collisions in the dilute limit. Few hydrodynamic calcu- 
lations are applied also to peripheral Pb-Pb collisions at 
= 2.76TeV [l,[2l,[2l| with resuhs compatible with ex- 
perimental observations. Nevertheless, it must be noted 
that as the impact parameter increases, uncertainties of 
the hydrodynamic model become more important; fiuc- 
tuations modify substantially the initial eccentricity, the 
relative role of the hadronic corona in the evolution of 
the system increases. In the present calculation the last 
issue is partly taken into account through an increase 
of the shear viscosity to entropy ratio at lower tempera- 
tures. For d-Au colhsions at ^/snn — 200GeV the hydro- 
dynamic model is expected to break down, as indicated 
by the absence of jet quenching. However, there are no 
published experimental results concerning directly the el- 
liptic flow in d-Au collisions or estimates from hydrody- 
namic models at RHIC energies. 

Below we present results from event by event viscous 
hydrodynamic simulations for p-Pb and d-Pb collisions at 
y^sJfN = 4.4 and S.llTeV respectively. We use Glauber 
Monte-Carlo model initial conditions for the hydrody- 
namic evolution. We calculate particle spectra, charged 
particle pseudorapidity distributions, elliptic and trian- 
gular flow coefficients as function of pseudorapidity and 
transverse momentum. 



system 


TeV 


CNN 

mb 


dN 
dVPS 

at r/ps = 


p-Pb 


4.4 
8.8 


66.4 
73.4 


50 ±5 
65 ±9 


d-Pb 


3.11 
6.22 


63 
69.8 


80 ±5 
95± 10 



TABLE I. NN cross section (third column) and expected 
density of charged particles at mid-rapidity in the NN cm. 
(fourth column) for different systems and energies in central 
collisions. 
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FIG. 1. (Color online) The distribution of participant nucle- 
ons at different impact parameters (boxes) and the average 
number of nucleons as function of the impact parameter (solid 
line) for p-Pb interactions. 



II. SIZE AND SHAPE OF THE INITIAL 
FIREBALL 

The number of particles produced in a p-Pb or d-Pb 
interaction can be estimated from Npart the number of 
participant (wounded) nucleons in collision. The Glauber 
Monte-Carlo model generates a distribution of events 
with different source sizes (number of participant nucle- 
ons) and different shapes (distribution of participant nu- 
cleons in the transverse plane) . The binary collision con- 
tribution is expected to be numerically small. Moreover 
the number of binary collisions is roughly Npart — 1(2). 
The presence of a term depending on the number of bi- 
nary collisions cannot be separated from the functional 
dependence on Npart- The number of participant nu- 
cleons in the Glauber model depends on the NN cross 
section. 

p-Pb interactions at the LHC are planed at the cm. 
energy in the NN system starting at ^/snn = 4.4rTeV. 
This corresponds to proton and Pb momenta of 3.5TeV 
and 208 x I.38TeV, attainable with the present magnetic 
field configurations in the accelerator For deuteron 
beams it gives the energy ^/sNN = S.IITeV. The max- 
imal available NN cm. energy at the LHC is 8.8 and 
6.22TeV for p-Pb and d-Pb interactions respectively. For 



collisions of beams with different energies per nucleon, 
the NN cm. reference frame is shifted in rapidity with 
respect to the laboratory frame. The shift is ysh = 0.46 
and 0.12 for p-Pb and d-Pb interactions. All the cal- 
culations in the hydrodynamic model are made in the 
NN cm. frame. For the final emitted particles a boost 
is made to the laboratory frame to get spectra around 
mid-rapidity or pseudorapidity distributions. 

The NN cross section at different energies can be 
obtained from an interpolation of values at 200GeV 
2.76TeV and 7TeV [25, 26] {aNN ^ 42, 62, and 71mb 
respectively) using a formula of the form (tnn oc a -f- 
bln{y^SNN) + c\n^ {^/snn) ■ The resulting NN cross sec- 
tions from the Table U are used in our Glauber model 
calculation. We take a Wood-Saxon profile for the Pb 
nuclear density 

pix,y,z) = — ^° z r , (2.1) 

1 + exp U^/x^ + - RA)/a] 

with po 0.17fm"^, Ra = 6.55fm and a = 0.45fm, 
and an excluded distance for nucleons of 0.4fm, for the 
deuteron we use the Hulthen distribution (23] . 

Events at a given impact parameter are generated us- 
ing the GLISSANDO code for the Glauber model 
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FIG. 2. (Color online) The probability distribution of par- 
ticipant nucleons in p-Pb interactions. The three centrality 
classes considered in the simulations are defined by cuts in 
the number of participant nucleons. 
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FIG. 3. (Color online) Same as Fig. [2] but for d-Pb interac- 
tions. 



The distribution of participant nucleons at different im- 
pact parameters is shown in Fig. [T] for p-Pb interactions 
at 4.4TeV. We notice that the number of participant nu- 
cleons fluctuates strongly at a fixed impact parameter. 
The number of participant nucleons can be significantly 
above the average value (solid line in Fig. [1]). Defin- 
ing the most central collisions as a interval in the impact 
parameter is incorrect. The few percents most central 
events in terms of the number of participant nucleons 
[Npart > 18) have a participant multiplicity larger than 
the average Npart at zero impact parameter. The pic- 
ture is very similar for d-Pb collisions. In the experiment 
the centrality classes are defined by the track multiplic- 
ity, closely correlated with the number of participants in 
the model. In heavy-ion collisions the number of partic- 
ipants is correlated with the impact parameter. In p-Pb 
or d-Pb interaction it is preferable to define the central- 
ity classes for events using directly cuts in Npart- Figs. M 
and [3] show the probability density for events of a given 
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FIG. 4. (Color online) Eccentricity (solid line) and triangu- 
larity (dashed line) in p-Pb interactions as function of the 
number of participant nucleons. 
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for the two systems considered. For p-Pb events, 



we use three centraHty classes defined as 18 < A^. 
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FIG. 5. (Color online) Same as Fig. |4]but for d-Pb interac- 
tions. 



11 < Npart < 17, and 8 < Npart < 10 corresponding to 
centrality bins 0-4%, 4-32% and 32-49% out of all the 
inelastic events {Npart ^ 2). The unusual numbers for 
the centrality percentiles are fixed by the discrete variable 
Npart- For the d-Pb interactions, we choose 27 < Npart, 
16 < Npart < 26, and 10 < Npart < 15 corresponding to 
centrahty bins - 5%, 5 - 30% and 30 - 50%. 

The charged particle density at central pseudorapid- 
ity can be estimated from the multiplicity observed at 
a similar energy and for a similar number of partici- 
pant nucleons measured in peripheral Pb-Pb collisions 
at the LHC Q, interpolating the measured values of 
dN /drjps I < Npart /2 > at centralities 60 - 70% and 
70 — 80% to the average number of participant nucleons 
{Npart) corresponding to the most central bins consid- 
ered in p-Pb and d-Pb collisions. The energy depen- 
dence of dN/rjps is s^'^^ for p-p and s°'^^ for nucleus- 
nucleus collisions [2^. We take s^'^^ to extrapolate from 
y^SNN = 2.76TeV. The estimated values of the charged 
particle density at midrapidity are quoted in Table HI the 
uncertainty comes from the uncertainty in the measure- 
ments [9|] and in the value of the exponent in the energy 
dependence. 

The azimuthally asymmetric collective flow is driven 
by the asymmetry of the initial fireball. The initial ec- 
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FIG. 6. (Color online) Contour plot s{x,y,Tj^^ = 0) of the 
initial entropy density in a d-Pb collision with Npart = 24. 



centricity in events with Npart participant nucleons 
(Etr^'cos(2(0,-V^2))) 



£2 



< EST' > 



(2.2) 



is calculated in each event with respect to the eccentricity 
angle maximizing £2 j the sum runs over all participant 
nucleons at positions r.^, 0^, (...) denotes averaging over 
events. In a similar way the triangularity is 

and is calculated with respect the triangularity axis ■03 
in each event [3,[l9']. In Figs. |4] and [5] we plot the eccen- 
tricity and the triangularity as function of the number of 
participant nucleons in the fireball. In proton induced 
interactions, the eccentricity and the triangularity of the 
source are similar and decrease for central collisions. It is 
different for d-Pd collisions, the eccentricity is larger than 
the triangularity and increases for central events. The 
eccentricity in d-Pb interactions is caused by the asym- 
metric configuration of the two nucleons in the deuteron. 
Configurations with a large separation of the deuteron 
proton and neutron in the transverse plane have a large 
eccentricity and usually lead to a large number of par- 
ticipant nucleons in the Pb nucleus. This effect causes 
the increase of the eccentricity for the most central colli- 
sions in Fig. [SJ We note that the eccentricity in Glauber 
models can be modified by correlation effects pvl . [29j 

We assume that the initial entropy density in the fire- 
ball is proportional to the number of participant nucle- 
ons. The density in the transverse plane x, y is the sum 
of contributions from participant nucleons at positions 
Xi, Hi from the Pb nucleus N-{x, y) and from the proton 
N^{x,y) (or from the proton and/or the neutron in the 
deuteron) 



N±{x, y) = So V Gxp 



(x - x.,)^ + {v- yif 



2al 



The contribution from each nucleon is a Gaussian of 
width Gyj = 0.4fm. The final results show some depen- 
dence on the chosen value of ■ Using a smaller width of 
0.3fm/c we notice an increase by ~ 10% of the integrated 
elliptic and triangular flow for p-Pb collisions. A simi- 
lar effect has been observed in Ref. [l^]. The parameter 
So is fixed to reproduce the final multiplicity after the 
hydrodynamic evolution. The distribution in space-time 
rapidity ri\\ is asymmetric 

s{x,y, 7711 ) = /_ (7711 )N.{x,y) + /+ (77,1 )N+ {x,y) , (2.5) 

the profiles f± {ri\\ ) are of the form 



1± 



Vbe 



fivw) 



(2.6) 



ybeam is the beam rapidity in the NN cm. frame. 

The asymmetric emission in the forward (backward) 
rapidity hemisphere from forward (backward) going nu- 
cleons can be observed in the distribution of charged par- 
ticles in d-Au collisions at RHIC The distribution 
of the form (j2.5p has been used as the initial condition 
for the hydrodynamic evolution in modeling Au-Au col- 
lisions at RHIC yielding a satisfactory description of the 
directed flow [Slj. The parameters of the longitudinal 
profile 



fiv\\) = exp 



2^2 



(2.7) 



(2.4) 



the plateau width 2?7o and the width of the Gaussian tails 
are adjusted as initial conditions for 3 -I- 1-D viscous 
hydrodynamic calculations to reproduce the charged par- 
ticle pseudorapidity distributions in Au-Au collisions at 
200GeV (7/0 = 1.5, cr,, = 1.4 [3) and Pb-Pb collisions 
at 2.76TeV [Hi (770 = 2.3, = 1.4). For the present 
calculation we take cr^ = 1.4 and 7^9 = 2.35, 2.4 for in- 
teractions at 3.11 and 4.4TeV respectively. An example 
of the initial entropy density in a d-Pb interaction event 
in shown in Fig. |6l Typically we observe strongly de- 
formed lumpy initial states. The elongated shape of the 
source results from the configuration of the nucleons in 
the deuteron while hitting the larger nucleus. This con- 
figuration is more important for the resulting eccentricity 
and the total number of participant nucleons than the im- 
pact parameter (as long as the deuteron hits the core of 
the Pb nucleus). 



III. VISCOUS HYDRODYNAMICS 

We use the second order relativistic viscous hydrody- 
namics to evolve the initial energy density in each event 
[3^ . The initial entropy density is generated in the 
Glauber Monte-Carlo procedure described in the previ- 
ous section. The viscous hydrodynamics incorporates de- 
viations from local equilibrium in terms of the shear and 
bulk viscosities, at zero baryon density heat conductivity 
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can be neglected. These corrections tt^'^ and 11 to the 
energy momentum tensor T^'^ are evolved dynamically 

AM"A'''^M^(9^7ra^ = -7^^"'^a,u" (3.1) 

Ttt 3 



and 



'Cd^u'' - n 4 
Tn 3 



1 



and = A'^'^di,. The hydrodynamic equations 



d^T^'" = 



(3.2) 



(3.3) 



(3.4) 



are solved numerically in the proper time t = ^/t^ — z'^ 
on a grid in the transverse coordinates x, y and the space- 
time rapidity ryu, starting from tq = 0.6fm/c. We use 
So = 0.72GeV3 in (EH) for both p-Pb and d-Pb collisions, 
which gives the expected final multiplicities. We take for 
the relaxation time = and assume tu = t-^- The 
initial fluid velocity is taken as the Bjorken flow, the 
initial stress corrections from shear viscosity correspond 
to the Navier-Stokes formula, while the initial bulk vis- 
cosity corrections are zero n(ro) = 0. The details of the 
solution in 2 -I- 1-D and 3 -I- 1-D are given in [l^, • 

The shear viscosity to entropy ratio in our calculation 
is not constant. It takes the value rj/s = 0.08 in the 
plasma phase and increases in the hadronic phase [Tg'] 



s 



fHG{T) + {\- fHGiT)) 



VQGP 



(3.5) 



with tihg/s = 0.5, rjQGp/s = 0.08, and fHG{T) = 
l/(exp((r-r//G)/Ar) -M), where Thg = 130MeV, 
AT = 30MeV. The bulk viscosity is nonzero in the 
hadronic phase 



^(T) 



C,HG 



IdT) 



(3.6) 



with Chg/s ^ 0.04 and /^(r) 
l/(exp((r-rj/Ar(;)-Hl), where T^; = 160MeV, 
AT^ = 4MeV. The equation of state is an interpolation 
of lattice QCD results at high temperatures [3j] and a 
hadron gas model equation of state at lower tempera- 
tures. In constructing the equation of state we follow 
the procedure of [34|]. The temperature dependence of 
the sound velocity has no soft point [3 ■ 

The hydrodynamic evolution stops at the freeze-out 
temperature of 135MeV. At the freeze-out hypersurface 
particle emission is done following the Cooper-Frve for- 
mula in the event generator THERMINATOR [35[, with 
viscous corrections to the equilibrium momentum distri- 
bution /o 



f — fo+ shear + Sfbulk 



(3.7) 
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FIG. 7. (Color online) Constant temperature hypersurface 
r(r, X = 0,y, r^ii = 0) in a p-Pb interaction for the freeze- 
out temperature Tf = 135MeV (dashed line) and for IGOMeV 
(solid line). 



We use quadratic corrections in momentum for the shear 
viscosity 



Sfshe 



fo (1 ± /o) 



1 



2T2(e+p) 



P'^P^TT^, (3.8) 



and asymptotically linear corrections for the bulk viscos- 
ity 

Sfbulk = Cbuikfo (1 ± fo) (c^sU^P^ - (^"P^)' ~ I n 



3w^Pm 



with 



(3.9) 



Cbulk 



(27r)3 E 



(3.10) 

the sum runs over all the hadron species. 

In Fig. [7]is shown the freeze-out hypersurface at yyn — 
for a p-Pb event with Nj 



part 



24. The dense source 



survives for 5fm/c with the lifetime of the deconfined 
phase of 3.5fm/c (T > 160 MeV, solid line contour in 
Fig. El). 



IV. RESULTS 

For each centrality class 50 hydrodynamic event are 
calculated. For each event several hundred THERMINA- 
TOR events are generated and analyzed together. This 
reduces non-flow effects, which in this case come from 
resonances decays. The numerical gird for the hydrody- 
namic evolution is set in the NN cm. frame. The mo- 
menta of the emitted particles are boosted by Ush = 0.46 
and 0.12 for p-Pb and d-Pb collisions to obtain spectra 
in the LHC laboratory frame. 

The distribution of charged particles in pseudorapidity 
is shown in Fig. [5] for the three centrality classes defined 
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FIG. 8. (Color online) Pseudorapidity distribution of charged 
particles in p-Pb interactions at ^/snn = 4.4TeV. The 
dashed, dashed-dotted and solid lines correspond to the three 
centrality classes defined by the number of participant nucle- 
ons. The distributions are shown in the laboratory frame for 
the LHC experiments. 
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FIG. 9. (Color online) Same as Fig. [8] but for d-Pb interac- 
tions at ^ysMN ~ S.llTeV. The long-dashed line shows the 
pseudorapidity density for the most central p-Pb collisions at 
4.4TeV. 
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FIG. 10. (Color online) Transverse momentum spectra of 
Tr"*" in p-Pb interactions at ^sjviv ~ 4.4TeV, for y = in 
the laboratory frame. The dashed, dashed-dotted and solid 
lines correspond to the three centrality classes defined by the 
number of participant nucleons. 
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FIG. 11. (Color online) Same as Fig. [10] but for K+ . 



nucleons in the Pb nucleus j3(j|. For d-Pb collisions in 
the centrality bin Npart > 27 the particle multiplicity is 
similar as in 60 — 70% centrality Pb-Pb interactions. This 
makes the applicability of the hydrodynamic model even 
more justified in that case. The particle multiplicity in 



in Sec. |TT1 The density of charged particles at midrapidity 
for centrality bins Npart ^ 18 and 11 < Npart < 17 in p- 
Pb is larger than observed in Pb-Pb collisions at 2.76TeV 
for centrality 70—80%. One can expect a similar degree of 
collective acceleration as in peripheral Pb-Pb events. The 
particle multiplicity in p-Pb interactions is of the same 
order as in p-p interactions with the highest multiplicity 
analyzed by the CMS collaboration @ . While the nature 
of the high multiplicity p-p events is still unclear, the 
multiplicity in a p-Pb or d-Pb collision is simply related 
to the source size and density. 

The charged particle densities in pseudorapidity for p- 
Pb (Fig. |S]) and d-Pb collisions (Fig. ^ are asymmetric, 
reflecting the predominant emission from the participant 
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FIG. 12. (Color onhne) Same as Fig. [TO]but for d-Pb inter- 
actions at y/SNN = S.llTeV. 



FIG. 13. (Color online) Same as Fig. [Hbut for K+ . 



p-Pb events with Npart ^ 18 is similar as for d-Pb events 
with 16 < Npart < 26. We find, that in all the cases the 
charged particle density at midrapidity is to within 5% 
proportional to the number of participant nucleons. 

The transverse momentum spectra for tt"*" and K'^ 
emitted in p-Pb collisions are hardened by the collec- 
tive transverse flow generated in the hydrodynamic ex- 
pansion (Figs. [TUl and [TT|) . For more central collisions 
the spectra are slightly flatter, as more transverse flow 
is generated. A similar picture appears for transverse 
momentum spectra in d-Pb collisions. The spectra from 
Npart > 27 and 16 < Npart < 26 centrality bins have a 
similar effective slope and are harder than the ones for the 
10 < Npart < 15 centrality class. It is interesting to ob- 
serve that the transverse momentum spectra are harder 
in p-Pb than in d-Pb collisions. The transverse size of 
the fireball in p-Pb interactions is smaller but its density 
is higher, this leads to a faster transverse expansion than 
for d-Pb interactions. 

Fluctuating initial densities (Fig. [6]) have a nonzero 
eccentricity and triangularity (Fig. |4]and[5]). The short 
hydrodynamical expansion stage in these systems is suffi- 
cient to generate noticeable elliptic and triangular flows. 
Fig. [14] shows the pseudorapidity dependence of the p± 
integrated elliptic V2 and triangular flow coefficients 
in p-Pb interactions. In the calculations, 500 to 1500 
THERMINATOR events are generated from each hyper- 
surface generated in a 3 -I- 1-D viscous hydrodynamic evo- 
lution. The event plane orientations and are found 
and the elliptic V2 and triangular flow coefficients are 
calculated in each event. The average over the hydrody- 
namic events gives U2{5'2} = (^2) and t'sjvE's} — (^3). 
The second cumula nt flo w coefficients include flow fluc- 
tuations w„{2} = ^ {Vn)- A moderate dependence of the 
elliptic and triangular flows on centrality is seen. In p-Pb 
collisions both the eccentricity and the triangularity de- 
formations of the initial shape are fluctuation dominated. 
We observe some reduction of the collective flow at for- 
ward and backward pseudorapidities. This reduction is 
due to an increase of dissipative effects and a shorter 
life-time of the source at nonzero space-time rapidities 
[36l [37I . The form of the pseudorapidity dependence of 
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FIG. 14. (Color online) The elliptic and triangular flow coeffi- 
cients of charged particles as function of pseudorapidity in the 
laboratory frame in p-Pb interactions, for centralities — 4% 
(top panel), 4 - 32% (middle panel), and 32 - 49% (bottom 
panel). The dashed-dotted and dotted lines represent the el- 
liptic flow coefficients V2{2} and V2{'^2}, while the dashed 
and solid lines represent the triangular flow coefficients i'3{2} 
and walk's}. 



the harmonic coefficients of the flow in Fig. [14] must be 
taken with caution, as 3-f 1-D viscous hydrodynamic cal- 
culations cannot reproduce it accurately [l^ [l6| . The az- 
imuthally asymmetric flow is different in d-Pb collisions 
(Fig. ITSl) . The elliptic flow is signiflcantly larger than 
in p-Pb interactions, reaching 0.097 for central collisions. 
We notice a strong centrality dependence, V2 increases 
significantly for central collisions. The initial eccentricity 
of the source in d-Pb colhsions is large (Fig. [5]). The el- 
liptic flow fluctuations (the difference between W2{2} and 
i'2{^2} [3^) are relatively less important for the deuteron 
than for the proton induced interactions. The triangular 
flow in d-Pb is similar as in p-Pb collisions, and does not 
vary strongly with the centrality. 

The hydrodynamic response translates the initial de- 
formation of the fireball into the azimuthal asymmetry 
of the final flow in an event by event basis. The elliptic 
flow coefficient follows closer the initial eccentricity than 
the triangular flow follows the initial triangularity, both 
in the magnitude and the orientation of the event plane. 
This observation is in agreement with other studies j39| . 
The hydrodynamic response v^jt^ is larger for central 
collisions, where dissipative effects that reduce the flow 
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FIG. 15. (Color online) Same as Fig. [14] but for d-Pb inter- 
actions, for centralities — 5% (top panel), 5 — 30% (middle 
panel), and 30 — 50% (bottom panel). 
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FIG. 17. (Color online) Same as Fig. [T6]but for the triangular 
flow. 
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FIG. 18. (Color online) Same as Fig. \M but d-Pb inter- 
actions. The dashed, dashed-dotted and solid lines corre- 
spond to the three centrality classes defined by the number 
of participant nucleons, Npart > 27, 26 > Npart > 16 and 
15 > Npart > 10 respectively. 



system (Fig. [T5| is large, it increases significantly for 
central collisions, where a saturation of the dependence 
on the transverse momentum appears around IGeV. The 
triangular flovi' in d-Pb interactions (Fig. [TO]) is similar 
in magnitude as in p-Pb collisions, and shows almost no 



FIG. 16. (Color online) The elliptic flow coefficient of charged 
particles as function of transverse momentum around y = 
in the laboratory frame for p-Pb interactions. The dashed, 
dashed-dotted and solid lines correspond to the three cen- 
trality classes defined by the number of participant nucleons, 
Npart > 18, 17 > Npart > H aud 10 > Npart > 8 respectively. 



asymmetry are smaller (40j . 

The elliptic and triangular flow coefficients show a hy- 
drodynamic behavior as function of the transverse mo- 
mentum (Figs. [TBI and ITS)) . Same as for the integrated 
flow, there is little change with centrality for the ellip- 
tic flow, while some decrease of in most central p-Pb 
collisions is seen. The elliptic flow V2{p±) for the d-Pb 




P^ [GeV] 

FIG. 19. (Color online) Same as Fig. I18l but for the triangular 
flow. 



9 



0.08 r '■ 
0.06 r , * ALICE Pb-PbH 

0.04^ A —A ^ hydro p-Pb J 
■ hydro d-Pb 

0.02 r -. 

20 40 60 80 100 120 140 

dN/dri 
'ps 

FIG. 20. (Color online) Elliptic flow coefficient t'2{2} as func- 
tion charge particle density at central pseudorapidity. Hydro- 
dynamic calculations for p-Pb collisions at y^sjvjv = 4.4TeV 
(triangles), for d-Pb collisions at S.llTeV (squares), and ex- 
perimental data from the ALICE collaboration for Pb-Pb col- 
lisions at 2.76TeV (stars) [JJ] are shown. 



variation with the centrality. 



V. CONCLUSIONS 

The formation of a hot, collectively expanding fireball 
in p-Pb collisions at ^snn = 4.4TeV and d-Pb collisions 
at S.llTeV is studied. We perform 3-1- 1-D event by event 
viscous hydrodynamic calculations. The initial size and 
shape of the fireball is taken from the Glauber Monte- 
Carlo model. The initial entropy density is adjusted to 
reproduce the expected particle multiplicity estimated as 
an extrapolation from observations in peripheral Pb-Pb 
collisions at ^/snn — 2.76TeV. 

A small hot and dense fireball is formed. It expands 
rapidly in the transverse direction. The p± spectra of 
emitted particles get harder, especially for p-Pb colli- 
sions. The deconfined phase survives for 3 — 4fm/c in 
events with high particle multiplicity, the presence of 
such a dense medium should be visible in the nuclear 
attenuation factor for high p± hadrons. The size and 
the life-time of the source can be further constrained in 
same-pion interferometry measurements. The initial ec- 



centricity and triangularity of a lumpy initial fireball lead 
to the formation of an azimuthally asymmetric flow. The 
elliptic flow is 3 — 4% in p-Pb collisions, with little cen- 
trality dependence (Fig. [^0]) . For the d-Pb system, the 
elliptic flow is significantly larger, increasing for central 
collisions, and reaching almost 10%. A comparison to 
peripheral Pb-Pb collisions in Fig. [SD] shows that simi- 
lar conditions are realized in proton or deuteron induced 
interactions. The elliptic flow of that magnitude can be 
measured, with a different dependence on centrality in 
p-Pb and d-Pb collisions. 

Let us close with a discussion on future prospects for 
proton and deuteron induced reactions at ultrarelativis- 
tic energies. p-Pb collisions at 4.4TeV are planned in the 
near future at the LHC Q. The elliptic flow and the 
hardening of the pj_ spectra are noticeable and should 
be looked for in the experimental analysis. However, it 
must be stressed that the dynamics of such small systems 
is at the limit of the applicability of the viscous hydro- 
dynamic model. The use of the hydrodynamic model 
in d-Pb interactions is better justified, also the elliptic 
flow is stronger, but such collisions are not planned in 
the near future at the LHC. The shift to the maximum 
LHC energy to = 6.22 for d-Pb and 8.8TeV for 

p-Pb collisions results in an increase of the particle mul- 
tiplicity by 30%. Hydrodynamic expansion would last 
longer, with less dissipative effects. The eccentricity and 
triangularity are similar as at lower LHC energies and 
qualitatively we expect similar results. 

In view of the results in this paper it seems very in- 
teresting to look for collective effects in d-Au collisions 
at ^/Jn = 200GeV in RHIC experiments. The multiplic- 
ity in central d-Au interactions is similar as in peripheral 
Au-Au collisions at the same energy. If some stage of col- 
lective expansion is present, the large initial eccentricity 
in a d-Au system should translate into a measurable ellip- 
tic flow. Unfortunately no published data exist for these 
experiments, hydrodynamical simulations are underway. 
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